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We consider a small metallic particle (quantum dot) where ferromagnetism arises as a consequence 
of Stoner instability. When the particle is connected to electrodes, exchange of electrons between 
the particle and the electrodes leads to a temperature- and bias-driven Brownian motion of the 
direction of the particle magnetization. Under certain conditions this Brownian motion is described 
by the stochastic Landau-Lifshitz-Gilbert equation. As an example of its application, we calculate 
the frequency- dependent magnetic susceptibility of the particle in a constant external magnetic field, 
which is relevant for ferromagnetic resonance measurements. 
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I. INTRODUCTION 

The description of fluctuations of the magnetization 
in small ferromagnetic particles pioneered by Browni is 
based on the Landau-Lifshitz-Gilbert (LLG) equation^ 
with a phenomenologically added stochastic term. This 
approach has been widely used: just a few recent appli- 
cations are a numerical study of the dynamic response 
of the magnetization to the oscillatory magnetic field;^ 
a numerical study of ferromagnetic resonance spectra^^ 
study of resistance noise in spin valves,^ and a study of 
the magnetization switching and relaxation in the pres- 
ence of anisotropy and a rotating magnetic fieldfi 

In equilibrium the statistics of stochastic term in 
the LLG equation can be simply written from the 
fluctuation-dissipation theorem^i However, out of equi- 
librium a proper microscopic derivation is required. Mi- 
croscopic derivations of the stochastic LLG equation out 
of equilibrium, available in the literature, use the model 
of a localized spin coupled to itinerant electronS )^i^i^°i^^ 
or deal with non-interacting electrons^ In contrast to 
this approach, we start from a purely electronic system 
where the magnetization arises as a consequence of the 
Stoner instability. Our derivation has certain similarity 
with that of Ref. 13 for a bulk ferromagnet, where the di- 
rection of magnetization is fixed and cannot be changed 
globally, so its local fluctuations are small and their de- 
scription by a gaussian action is sufficient. This situation 
should be contrasted to the case of a nanoparticle where 
the direction of the magnetization can be completely ran- 
domized by the fluctuations, so that the effective action 
for the direction of the classical magnetization has a non- 
gaussian part. The bias-driven Brownian motion of the 
magnetization with a fixed direction (due to and easy- 
axis anisotropy and ferromagnetic electrodes) has been 
also studied in Ref. [13 using rate equations. 

We assume that the single-electron spectrum of the 
particle, which is also called a quantum dot in the litera- 
ture, to be chaotic and described by the random-matrix 
theoryi^ii^. To take into account the electron-electron in- 



teractions in the dot we use the universal Hamiltoniaufii 
with a generalized spin part, corresponding to a ferro- 
magnetic particle. Electrons occupy the quantum states 
of the full Hamiltonian and form a net spin of the par- 
ticle of order of 5*0 ^ 1; throughout the paper we use 
h — 1. The dot is coupled to two leads, see Fig. (TJ which 
we assumed to be non-magnetic. The approach can be 
easily extended to the case of magnetic leads. The num- 
ber Nch of the transverse channels in the leads, which are 
weU coupled to the dot, is assumed to be large, Nch ^ 1- 
Equivalently, the escape rate 1/r of electrons from the 
dot into the leads is large compared to the single-electron 
mean level spacing Si in the dot. This coupling to the 
leads is responsible for tunnelling processes of electrons 
between states in the leads and in the dot with random 
spin orientation. As a result of such tunnelling events, 
the net spin of the particle changes. We show that this 
exchange of electrons gives rise both to the Gilbert damp- 
ing and the magnetization noise in the presented model, 
and under conditions specified below, the time evolution 
of the particle spin is described by the stochastic LLG 
equation. 

We study in detail the conditions for applicability of 
the stochastic approach. We find that these limits are 
set by three independent criteria. First, the contact re- 
sistance should be low compared to the resistance quan- 
tum, which is equivalent to Nch S> 1. If this condi- 
tion is broken, the statistics of the noise cannot be con- 
sidered gaussian. Physically, this condition means that 
each channel can be viewed as an independent source 
of noise, so the contribution of many channels results in 
the gaussian noise by virtue of the central limit theo- 
rem if Nch S> 1. Second, the system should not be too 
close to the Stoner instability: the mean-field value of 
the total spin 5*0 ^ \/ Net- If this condition is violated, 
the fluctuations of the absolute value of the magneti- 
zation become of the order of the magnetization itself. 
Third, 5*0 ^ T^s/Si, where T^e is the effective tempera- 
ture of the system, which is the energy scale of the elec- 
tronic distribution function determined by a combination 
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FIG. 1: (Color online). Device setup considered in this work: 
a small ferromagnetic particle (quantum dot) coupled to two 
non-magnetic leads (see text for details). 

of temperature and bias voltage (the Boltzmann constant 
/cb = 1 throughout the paper). Otherwise, the separation 
of the degrees of freedom into slow (the direction of the 
magnetization) and fast (the electron dynamics and the 
fluctuations of the absolute value of the magnetization) 
is not possible. 

In the present model we completely neglect the spin- 
orbit interaction inside the particle, whose effect is as- 
sumed to be weak as compared to the effect of the leads.— 
The effects of the electron-electron interaction in the 
charge channel (weak Coulomb blockade) are suppressed 
for Nch ^ Ifi^ so we do not consider it. 

As an application of the formalism, we consider the 



magnetic susceptibility in the ferromagnetic resonance 
measurements, which is a standard characteristic of mag- 
netic samples. Recently, a progress was reported in mea- 
surements of the magnetic susceptibility on small spatial 
scales in response to high-frequency magnetic fields.^— 
Measurements of the ferromagnetic resonance were also 
reported for nanoparticles, connected to leads for a some- 
what different setup in Ref. [20l . 

The paper is organized as follows. In Sec. |TT]we intro- 
duce the model for electrons in a small metallic particle 
subject to Stoner instability. In Sec. IIIII we analyze the 
effective bosonic action for the magnetization of the par- 
ticle. In Sec. lIVI we obtain the equation of motion for the 
magnetization with the stochastic Langevin term, which 
has the form of the stochastic Landau-Lifshitz-Gilbert 
equation, and derive the associated Fokker-Planck equa- 
tion. In Sec.|V]we discuss the conditions for the applica- 
bility of the approach. In Sec. IVII we calculate the mag- 
netic susceptibility from the stochastic LLG equation. 



II. MODEL AND BASIC FORMALISM 

Within the random matrix theory framework, elec- 
trons in a closed chaotic quantum dot are described by 
the following fermionic action: 



dt 



N 



n,n' — l 



N 

«=1 



(2.1) 



Here ipn is a two-component Grassmann spinor, t runs along the Keldysh contour, as marked by f ; a^-y'^ are the Pauli 
matrices (we use the hat to indicate matrices in the spin space and use the notation cto for the 2x2 unit matrix). 
Hnn' is an N X N random matrix from a gaussian orthogonal ensemble, described by the pair correlators: 



N5l 



7nn' ^7im' 



(2.2) 



Here 5i is the mean single-particle level spacing in the dot. 

The magnetization energy E{S) is the generalization of the JsS'^ term in the universal Hamiltonian for the electron- 
electron interaction in a chaotic quantum dotJ^ Since we are going to describe a ferromagnetic state with a large value 
of the total spin on the dot, we must go beyond the quadratic term; in fact, all terms should be included. E{S) can be 
viewed as the sum of all irreducible many-particle vertices in the spin channel, obtained after integrating out degrees 
of freedom with high energies (above Thouless energy); the corresponding term in the action is thus local in time, and 
can be written as the time integral of an instantaneous function E{S{t)). This functional can be decoupled using the 
Hubbard- Stratonovich transformation with a real vector field h(t), which we call below the internal magnetic field: 



exp (^i j> dtE{S)^ = J Vh{t) exp j> dt{2h ■ S - Eih)) 
We rewrite the action in the form 



Slip.ip^h] = idt 



N 



V't(t)(G-M ^i:,At)-E{h{t)) 



(2.3) 



(2.4) 
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where the inverse Green's function 

( G^M = {i&odt + h-&) 6nn' - Hnn'^O (2.5) 
\ / nn' 

is a matrix in time variables t, t' , in orbital indices n and 
n' with 1 < n, n' < iV, in spin indices, and in forward (+) 
and backward (— ) directions on the Keldysh contour. 
Integration over fcrmionic fields ipn , V'i yields the purely 
bosonic action: 

S[h] = -iTr \ln{-iG-^)^ - dt E{h{t)), (2.6) 

where the trace is taken over all indices of the Green's 
function, listed above. 

In the space of forward and backward directions on 
the Keldysh contour, we perform the standard Keldysh 
rotation, introducing the retarded (G^), advanced {G^), 
Keldysh (G^), and zero (G^) components of the Green's 
function: 



1/1 1 Vg++ g+- 

2 I 1 -1 [ G-+ G- 



1 1 

-1 1 



( G^ G^ 
\G^ G^ 

(2.7) 

as well as the classical (/i^') and quantum {h'^) compo- 
nents of the field: 



h'^i n't \ If 1 -l \ ( h+ 
h^t /i^W ^ 2 I 1 1 M h 



1 1 

1 -1 



(2.8) 

We will also write this matrix as /i = /I'^'r^' + /I'^r'?, 
where r^' and are 2x2 matrices in the Keldysh space 
coinciding with the unit 2x2 matrix and the first Pauli 
matrix, cTx, respectively. 

The saddle point of the bosonic action Eq. (|2.6p is 
found by the first order variation with respect to h'^''''^{t), 
which gives the self-consistency equation: 



hi(t) = 0, 



dEih^^it)) _ i 
dhf{t) ~ 2 



Tr{a^Gl(t,t)}. (2.9) 



We also note that the right-hand side of this equation is 
proportional to the total spin of electrons of the particle 
for a given trajectory of /i'^'(i): 

S(t) = ^Tr{a^G^(^,i)}. (2.10) 



In Eqs. (|2.9p and (|2.10p . the trace is taken over orbital 
and spin indices only. 

In the limit iV — s- cx), one can obtain a closed equation 
for the Green's function traced over the orbital indices;-^ 



N 



g(t,t') = ^y^Gnnit,t'). (2.11) 

n=l 

The matrix g{t,t') satisfies the following constraint: 

/ git, t") git", t') dt" = T^'a^5{t - t'), (2.12) 



where the right-hand side is just the direct product of 
unit matrices in the spin, Keldysh, and time indices. 
The Wigner transform of g^it,t') is related to the spin- 
dependent distribution function fie,t) of electrons in the 
dot: 



g^it + i/2, t - i/2) e*^' di = 2}{e, t) 



(2.13) 



In equilibrium, /(e) = ctq tanh(£/2T). 

The self-consistency condition (|2.9p takes the form 



dEih^^t)) _ TT 



— limTr{a^g^(M')} 



2hfit) 



dhfit) 26i i'^'t r ^" ^ Si 

(2.14) 

The last term takes care of the anomaly arising from 
non-commutativity of the limits N oo and t' — > t. 

In this paper we consider the dot coupled to two leads, 
identified as left (L) and right (i?). The leads have 
Nl and Nr transverse channels, respectively, see Fig[T] 
For non-magnetic leads and spin-independent coupling 
between the leads and the particle, we can characterize 
each channel by its transmission T„ with < T„ < 1 
and by the distribution function of electrons in the chan- 
nel Tnit — t'), assumed to be stationary. We consider 
the limit of strong coupling between the leads and the 
particle, Y.n=iTn » 1. 

The coupling to the leads gives rise to a self-energy 
term, which should be included in the definition of the 
Green's function, Eq. (|2.5p . Without going into details 
of the derivation, presented in Ref. [U, we give the final 
form of the equation for the Green's function traced over 
the orbital states, Eq. (|2.1ip : 



[dt - ih ■ a-,g] 



^ 27r 

n— 1 



V- TJi f -T,,g^ gRTn - T^g^ - g 



Tn 

2 



n 
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K 



+ Tng^ g^Tn + Tng^ 
-ff-^ - i + g^:Fn 



(2.15) 



Here the products of functions include convolution in 
time variables. This equation is analogous to the Usadel 
equation used in the theory of dirty superconductors 

To conclude this section, we discuss the dependence 
E{h). Deep in the ferromagnetic state, i.e. far from the 
Stoner critical point, we expect the mean-field approach 
to give a good approximation for the total spin of the 
dot. Namely, the mean field acting on the electron spins, 
is given by 2/io = dEiS)/dS = E'iS). We then require 
that the response of the system to this field gives the 
same average value for the spin: 



5n - 



2/io E'iSo) 



25i 



25i 



(2.16) 



Here we evaluated 5*0 from Eq. (|2.10p and applied the 
self-consistency equation (j2.14p to equilibrium state with 
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oc i7o, when the contribution of the first term in the 
right hand side of Eq. (|2.14p vanishes. 

Not expecting strong deviations of the magnitude of 
the spin from the mean-field value, we focus on the form 
of E{h) when \h\ k, Hq. The inverse Fourier transform of 
Eq. (|2.3p and angular integration for the isotropic E{S) 
gives 



^-^E(k)At ^ ^^^^^ I ^^^J^ShAt ^-^EiS)Ats2^s^ (2.17) 



where Ai is the infinitesimal time increment used in the 
construction of the functional integral in Eq. (|2.3D . 
Expanding E{S) near the mean- field value 5*0, 

E{S) « E{Sq) + E'{S - 5o) + ^ (5 - ^o)', (2.18) 

performing the integration in the stationary phase ap- 
proximation and using Sq = Kq/ 5i — — £"/(2(5i), we ob- 
tain 

where E^ is /i— independent term. This expression for 
E{h) defines the action S[h], Eq. (j^ . 

The energy £'(5') does not contain the energy i?oib(5'), 
associated with the orbital motion of electrons in the 
particle. Namely, to form a total spin S of the par- 
ticle, we have to redistribute 5 electrons over orbital 
states, which changes the orbital energy of electrons by 
i?orb(5') ^ SiS'^. The total energy £'tot(S') of the particle 
is the sum of two terms: Etot{S) — E{S) + Eorh{S). Sim- 
ilarly, we obtain the total energy of the system in terms 
of internal magnetic field 



Etot{h)^E{h) 



Si 

1 1 



(2.20) 



[h-haf+Ei, 



where Ei does not depend on h. We notice that the 
extremum of Etot (h) corresponds to h = ho and describes 
the expectation value of the internal magnetic field in an 
isolated particle. The energy cost of fluctuations of the 
magnitude of the internal magnetic field is characterized 
by the coefHcient 1/E" + l/26i. 



III. KELDYSH ACTION 

In this Section we analyze the action Eq. (|2.6p for the 
internal magnetic field h". We expect that the classical 
component ft.'^' (t) of this field contains fast and small os- 
cillations of its magnitude around the mean-field value 
ho. We further expect that the orientation of h'^\t) 
changes slowly in time, but is not restricted to small 



deviations from some specific direction. Based on this 
picture, we introduce a unit vector n{t), assumed to de- 
pend slowly on time, and write 



h^'it) = iho + hf{t))n{t), 



(3.1) 



where ft.|j'(i) is assumed to be fast and small. We expand 
the action (|2.6p to the second order in small fluctuations 
of the quantum component h"^ (t) and the radial classical 
component /i|j'(i): 

S[h] ~ -— I dtg''{t,t)h'>{t) + 

+ ^ j dthl\t)nit)h^t) 

dtdt' Ufj{t,t')hj{t) hl'{t')nj{t') 



- / dtdt' nf^{t,t') hf{t)n,{t) h'^^it') 



dtdt'nfjt,t')hfit) M{t' 



(3.2) 



The applicability of this quadratic expansion is discussed 
in Sec. IVBl 

In Eq. (|3.2p we introduced the polarization operator, 
defined as the kernel of the quadratic part of the action 
of the fluctuating bosonic fields: 



u^i^lt t') = 



Yiq,cl Yiq,q 

i J^TrjlnG-i} 



2 5h'^{t')5hf{t) 



(3.3a) 



(3.3b) 



where a,(3 — cl,q and i,j = x,y,z. The short time 
anomaly is explicitly taken into account in the definition 
of the polarization operators, see Eq. (|3.14p below. 

The first term of Eq. (|3.2p contains the vector of the 
Keldysh component of the Green function g 



K 



K 



K 



(T ■ g'^ . We emphasize that the Green's function and 
the polarization operator in Eq. (|3.2p are calculated at 
/ijj'(i) — and h'^{t) = for a given trajectory of the 
classical field hon{t). 

A. Keldysh component of the Green function 

For the Green's function in the classical field we have 
g''{tX) = -9^{t.t')^a%{t^t'), (3.4) 
while the Keldysh component satisfies the equation 

9")- (3.5) 



[^,-^hon■a,g^]^J2^i'^^" 



n=l 

We introduce the notation 



l^gT„(5i 1 , 1 



ri,=l 



27r tl tr 



(3.6) 
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Then the scalar and vector components of 



K 



& ■ satisfy two coupled equations: 



dt' + - 

T 



^ — ' ZTT 

n=l 

+tho[n{t)-n{t')]-g''{t,t'), 
1 



(3.7a) 



dt + df + 



g''{t,t') + ho[n{t)+n{t')]xg''{t,t') 



= tho[n{t)-n{t')]g^{t,t') 



(3.7b) 



As a zero approximation, we can consider the station- 
ary situation: gQ{t,t') — g^{t — t') and n{t) = const. In 
this case, we have 



9oM 



TL 



■2TL{t't') + —2TR(t~t') (3.8) 



and g^ — 0. 

For an arbitrary time dependence n(t), Eq. (I3.7bp can- 
not be solved analytically. However, if the variation of 
n{t) is slow enough, we can make a gradient expansion: 



,K 



+ 2hon x g-*^ = ithohg^ . (3.9) 



Here we introduced t = {t + t')/2, i = t — t' , dt + dt> = di- 
The dependence on t is split off and remains unchanged, 
while for the dependence on t the solution is determined 
by a linear operator £+: 



1 



±di± 2han x 



(3.10a) 



C+{u;)X{u;) 
1 



n{n ■ X{u!)) 
—iuj + 1/r 
-n X [n X X{uj)] ± i[n x X(uj) 



'i{u;±2ho) + l/i 



(3.10b) 



Here we assume that the direction of the internal mag- 
netic field n changes slowly in time, and \ri\T <C 1. 

Thus, all perturbations of g^ decay with the charac- 
teristic time T. In particular, the solution of Eq. p.9p 
has the form 

K ~ + K ~ f^- n + 2hoT n X h 
g =ithoC:^ngQ k. ithoT g^ ^ ^ — . (3.11) 



(2/ior)2 



1 



Expression for the first term in Eq. p.2p can be easily 
obtained from Eq. p. lip by taking the limit t — > and 
taking into account that any fermionic distribution func- 
tion in the time representation has the following equal- 
time asymptote: 



1 



int-t' 



{t^t'). 



We have 



2hoT n + 2hQTn x n 



(2/ior) 



1 



(3.12) 



(3.13) 



We notice that n ■ g^{t,t) = 0, and therefore the first 
term in the action Eq. (|3.2p is coupled only to the tan- 
gential fluctuations of h'^{t) _L n{t). 



B. Polarization operator 

We express the polarization operator in terms of the 
unit vector n{t). The polarization operator can be rep- 
resented as the response of the Green's functions to a 
change in the field, as follows directly from the defini- 
tion (j3.3b|l and the expression (|2.6p for the action: 



(3.14) 

Here the Green function Sg{t, t) can be calculated as the 
first-order response of the solution of Eq. I|2.15p to small 
arbitrary (in all three directions) increments of 5h'^\t) 
and 5hf{t). The zero-order solution of Eq. (|2.15p in the 
field h'^^ = hoTi and h,^ = is 



First, we calculate 6g^ , which responds only to dh"^: 



dt + dt' 



5g^{t, t') - iho& ■ n{t) dg^{t, t') 



+6g^{t, t') iho& ■ n{t') = 2icr ■ 5hi{t) 5{t - t'). 

(3.16) 

Since dt + df — c%, the solution always remains propor- 
tional to (5(t — t'): 



5f{t, t') = -2ia{C-5h'i){t) 5{t - t'). 



(3.17) 



Given 5g^ , components 6g^'^ can be found either from 
Eq. (|2.15p . or, equivalently, using the constraint g^ = 1: 
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Z -K 



g6g + 6gg^0 ^ 6g^ ^ 2 ' 2 ' 

(3.18) 

We notice that both Sg^'^ respond only to h'^{t) and, 
therefore, 



TT{aHSg"it,t)+dg^{t,t))} 
5hf{t') 



0. (3.19) 



This equation ensures that the action along the Keldysh 
contour vanishes for /i' = 0. 

To evaluate the remaining three components of the po- 
larization operator, we can apply the variational deriva- 
tives to the sum of 6g^{t,t) + 6g^{t,t) with respect to 
either classical 5hf'^{t') or quantum 5hfl{t') field, which 
give Hg(t,f') and Wfj{t,t'), respectively. Then, the ad- 
vanced component H^(t,t') = [\lf^{t' ,t)]* . 
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The equation for 5g^ = cr ■ Sg^ reads as 
1' 



dt + dt> + 



+ ho [n{t) X 5g^{t, t') - Sg'^it, t') x n{t')\ 
i[6h'\t)~5h'\t')\g^{t^t') 
-2i6h''{t)6{t-t')-Q{t,t'), 

(3.20a) 



where 
Q{t,t') 



n=l 

E 



27r V 2 2 

T„(l - Tn)Sl 



2tt 



9o 



- 5g 



90 



(3.20b) 



and ig^ = Tr{o-(55^}/2 with bg^ given by Eq. (|3T7| . 

To calculate the retarded component of the polar- 
ization operator, we calculate the response of Sg^{t,t') 
to (5/i^ in the limit t' — + t. Using the asymptotic behavior 
of the Fermi function, Eq. (|3.12p . we obtain: 



(3.21) 



Substituting this expression for 6g^{t,t) to Eq. (|3.14p . 
we obtain 



n^(^) 



n 



(co)+nf ,,(^), (3.22a) 



with 



n5,,(^) 



Z niTij 
2 v-^ ~ 



(3.22b) 
(3.22c) 



{l±2ihQT) 

1 — {{uj =F 2hQ)T ' 

Here we represented the polarization operator 11^ (tj) as 
a sum of the radial, nj|^-^-(a;), and tangential, U^^^jluj), 
terms. We note that the action Eq. (|3.2p contains only 
the radial component of the retarded and advanced polar- 
ization operators because we do not perform expansion 
in terms of the tangential fluctuations of the classical 
component of the field h,'^' (t) . 

In response to Sh"^, both corrections dg^{t,t') and 
5g^{t,t') contain terms oc d{t — t'), However, their sum 
6g^ {t, t') + 5g^ {t, t') remains finite in the limit t t': 



Sg''{t,t') + Sg^{t,t')^~2i / dt" / dh dt2 C+{t - h) Q{ti - t2 + i/2;t2 - h + i/2) C-{t2 ~ t") Sh''{t"), (3.23a) 



Q{T,;r2)^-S{n)6{r2)~J2^"^' 



T 

E 

n=l 



71=1 



27r 



+ ^n(Tl)^n(T2) 



(3.23b) 



r„(i-T„)^i 

27r 



5(f(n) 



9^ir2) 



- ^n{T2) 



with t^{t + t') /2 and i ^ t - t' . 



imation of non-interacting electrons 



From Eq. (|3.23|) we obtain the following expression for 
the Keldysh component of the polarization operator: 



71 — 1 

X { [8~g^{e)g^{e + Lo)~ATr.{e)^n{e + Lu)] 

+ (1 - T„) [gl^ie) - 2Tn{e)] [gl^ {e + uj) - 2Tn{e + u;)] }. 

(3.25) 

In principle, electron-electron interaction in the charge 
channel can be taken into account. The interaction mod- 
ifies the expression Eq. (|3.25p for TZ{uj) to the higher 
order^'^ in tSi ^ 1 and we neglect this correction here. 

In this paper we consider a particle connected to 
electron leads at temperature T with the applied bias 
V. In this case, Tl,r{£) = tanh(£ — /zl,_r)/(2T) with 



n,^(c.) = Hf,^.(c.) + nf,,,^.(c.), 



rLiTij 



w2 + l/r- 



(3.24a) 
(3.24b) 



Here function 72, (w) coincides with the noise power of 
electric current through a metallic particle in the approx- 



7 



^■L — I^R = V, and the integration over e gives 



and V^iuj) = [D^{lj)]*. The Keldysh component is 



27rT7^(w) =4LJCoth— + STt(F,w), (3.26) 



where 



Tt{V, w) = ^ 2(u; ± l^) coth - 4w coth ^ 



2T 



and ^ is the "Fano factor" for a dot 



^ nS-R 



rn(i - T„). 



2T 
(3.27) 



(3.28) 



At 3> T the function Tr(F, w) has two scales of u: 
(i) T smears the non-analyticity at (jJ — > 0, but the value 
of Tt{V,lo) deviates from Tt[V,Q) at \lu\ - \V\. Thus, 
the typical time scale above which one can approximate 
li^ {ijj) by a constant is at least lo <C max{r, \V\]. In the 
limit cij ^ we have 

n, (c = 0) = ^ (^n.n, + j^j^^;^^ ) ■ (3-29) 

The effective temperature Toff is given by 

Toff ^ T + coth ^-tV (3.30) 



C. Final form of the action 

We can rewrite the action for magnetization field h = 
{h'^'';h'^} with h'^' in the form of Eq. (|3.ip as a sum of 
the radial and tangential terms: 

S[h] ^Si\[hf,h1^]+S±[n{t),hl]. (3.31) 

The radial term in the action has the form 

5|| [hf, hi] = [V-'r^it, t') h^{t) /ij(i'), (3.32) 

where the inverse function of the internal magnetic field 
propagator is given by 



Tlf{t,t') 



(3.33) 



and Il"^^{t,t') — riinjll"^^. (t, t'). From this equation we 



find 



I?f(^)=I?f'(c.) 



-iuj + 1/r 



4 ~iuj+{5i+E"/2)/{T5^y 

(3.34) 



Vf{uj) = V^^'^iu) = Vf{uj)I{f{u)V^{uj). (3.35) 
The tangential term in the action is 



^1 J - (2/ior)2 + l 
^ [ dt [n{t) X [n{t) X B]] ■ h\{t) 

01 J 

dtdt'hl/t)ul,^it~t') hl^it'). 

(3.36) 



Here we recovered the external magnetic field B{t). The 



polarization operator 11^ j,- is given by Eq. p.24cp . 



IV. LANGEVIN EQUATION 

A. Langevin equation for the direction of the 
internal magnetic field 

In this section we consider evolution of the direction 
vector n, described by the tangential terms in the action, 
Eq. (|3.36p . We neglect fluctuations of the magnitude 
of the internal magnetic field, /ly, the conditions when 
these fluctuations can be neglected are listed in the next 
section. 

We decouple the quadratic in h.^ component of the 
action in Eq. (I3.36P by introducing an auxiliary field 'w{t) 
with the probability distribution 

V[w{t)] CX exp I ^ y" dtdt' {n^)T^\t,t') W^{t) Wj{t')^ , 

(4.1) 

and the correlation function 

{w.{t)w,{t'))^f^Ill,Jt,t'). (4.2) 



The field w{t) plays the role of the gaussian random 
Langevin force. Integration of the tangential part of 
the action, Eq. p.36p . over h"^^ produces a functional 
(5-function, whose argument determines the equation of 
motion: 

n + 2Tho\n X h] 1 , r , ^ 

4r4+l "x["x^])=0. (4.3) 

The above equation can be resolved with respect to h: 

h = -2[n X {w + B)] —[nx[nx{w + B)]] . (4.4) 

hoT 

This equation is the Langevin equation for the direc- 
tion n{t) of the internal magnetic field in the presence 
of the external magnetic field B{t) and the Langevin 
stochastic forces w{t). 
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B. The Fokker-Plank equation 

Next, we follow the standard procedure of derivation 
of the Fokker-Plank equation for the distribution V{n) 
of the probability for the internal magnetic field to point 
in the direction n. The probability distribution satisfies 
the continuity equation: 



dV dJi 



= 0, 



dt dn. 

where the probability current is defined as 

1 



(4.5) 



J 



2nx B 



hoT 



n X [n X B] rp 



(4.6) 



and the stochastic velocity ^ is introduced in terms of 
the field w as 

I = ~2[n X w]- [n X [n X w]] . (4.7) 

hoT 

The derivative d/dn is understood as the differentiation 
with respect to local Euclidean coordinates in the tangent 
space. Performing averaging over fluctuations of w in 
Eq. (|4.6p . we obtain 



dV _ d \ {2hoT)[n x B] + [nx[nx B] 
dt dn 



hoT 



1 d'^V 



(4.8) 



where the time constant Tq is defined as 



% = 



(4.9) 



Below we use the polar coordinates for the 
direction of the internal magnetic field, n = 
{sin0cosiy9,sin6'sini^,cos^?}. In this case the Fokker- 
Plank equation can be rewritten in the form 



dV 



I 



d 



dt sni9 dtp 
1 



Tq sin Q dip 



sm ( 



sin 6 FaV 



sine dV 



% de 



(4.10) 



where 



Bx sin Lp — By cos p 



hor 

+ 2cos9{Bx cos ip + By sin (p) - 2 sin 6* B^, (4. If) 
Fg = 2{Bx sin p — By cos p) 



cos 9 sin 9 
- (Bx cos p + By sm p) + _Bj, 



hoT 



hoT 



(4.12) 



It should be supplemented by the normalization condi- 
tion: 



2n TT 



j dp j sin9d9V{p,9) = 1, 



(4.13) 







which is preserved if the boundary conditions at = 0, tt 
are imposed: 



dV 

lim sin 9 I dp —— = . 



(4.14) 



Below we apply the Fokker Plank equation for calcu- 
lations of the magnetization of a particle 



M = I ^nV{n) 

Ait 



(4.15) 



V. APPLICABILITY OF THE APPROACH 

In this section we discuss the conditions of validity 
of the stochastic LLG equation, see Eq. (|4.f Op . for the 
model of ferromagnetic metallic particle connected to 
leads at finite bias. We briefly listed these conditions 
in the Introduction. Here we present their more detailed 
quantitative analysis. 



A. Fluctuations of the radial component of the 
internal magnetic field 

We represented the classical component of the internal 
magnetic field h'^^ in terms of a slowly varying direction 
n{t) and fast oscillations /ijj' of its magnitude around 
the average value /iq. Now, we evaluate the amplitude of 
oscillations of the radial component /ijj' of the field, using 
the radial term in the action, see Eqs. (|3.3f [) and (|3.32p . 

The typical frequencies for time evolution of small fluc- 
tuations of the internal magnetic field in the radial direc- 
tion are of order of 



5i + E"/2 I 



(5.1) 



as one can conclude from the explicit form of the prop- 
agator I?|j*(cL'), Eq. (|3.34p . of these fluctuations. This 
scale has the meaning of the inverse i?C-time in the 
spin channel. Deep in the ferromagnetic state (i. e., far 
from the Stoner critical point E" + 25i ~ 0) we estimate 
5i + E" 12 ~ bi (which is equivalent to E" ~ ho/So), 
so this spin-channel i?C-time is of the same order as the 
escape time r. This estimate for the frequency range is 
consistent with the simple picture, which describes the 
evolution of the internal magnetic field of the grain as 
a response to a changing value of the total spin of the 
particle due to random processes of electron exchange 



9 



between the dot and the leads. The electron exchange 
happens with the characteristic rate 1/r. 

The correlation function {hf{t)hf{t')) can be evalu- 
ated by performing the Gaussian integration with the 
quadratic action in and Using Eq. (|3.35p . we ob- 
tain the equal-time correlation function 



ll\2 



27r7^(w) 



(5.2) 



32rJi J 27r + [1 + E" /{25i)Y /t"^ 



This equation gives the value of fluctuations of the radial 
component of the internal magnetic field of the particle. 
These fluctuations survive even in the limit T = and 
y = 0, when 7^(ti>) — 2|w|/7rr. We have the following 
estimate 



{{h 



{E" 



■In 



Ett 



16ttt5i l + E"/{25i) 



(5.3) 



the upper cutoff E^ is the Thouless energy, E^ = vp/L 
for a ballistic dot with diameter L and electron Fermi 
velocity vp- 

The separation of the internal magnetic field into the 
radial and tangential components is justified, provided 

that the fluctuations 



; / of the radial component 

are much smaller than the average value of the fleld /lo, 
i.e. ((/i|jO^) < K- Using the estimate Eq. (jO)) . we 
obtain the necessary requirement for the applicability of 
equations for the slow evolution of the vector of the in- 
ternal magnetic field of a particle: 



5*0 > 



-^HE^t), 

TO I 



(5.4) 



where Sq is the spin of a particle in equilibrium and 
we again used the estimate E" ~ ho/So. Condition of 
Eq. (15. 4|) requires that the system is not close to the 
Stoner instability. 



At the same time, the typical value of /i|j'(a; ^ 1/t), 
as determined by the gaussian part of the action, was 
estimated in the previous subsection to be of the order of 

^ ^/t6i <C 1, so the higher-order terms 



are indeed not important. 

For the quantum component of the field the quadratic 
and quartic terms in the action are estimated as 

Sn+1 : / — X 



Si J {2t:Y 
X /i«(wi).../i?(w„)/i«(-cji 



■w„). (5.6) 



If Tcff 3> l/r, then the typical frequency scale is w ~ 1/t, 
so the quadratic term gives /i'^(u; ~ l/r) ^ y^fi/T^, and 

Sn - (<5l/reff)"/2-l ~ Hi/(TTeff)]"/2-l. If Teff « 1/t, 

at the typical scale u T^s we obtain h'^{uj ~ Tcff) ~ 
y/Si/iT^gT), SO again Sn ~ (t^i)"/2-i < 1 forn > 2. 

Physically, the parameter 1/{tSi) = Nch (or Teff/(5i, 
if it is larger) can be identified with the number of the 
independent sources of the noise acting on the magne- 
tization field. Thus, the smallness of the non-gaussian 
part of the action is nothing but the manifestation of the 
central limit theorem. 



C. Applicability of the Fokker-Plank equation 

From the above analysis we found that evolution of 
the direction of the internal magnetic field in time is 
described by a characteristic time Tq, introduced in 
Eq. (|4.9p . From the analysis of the fluctuations of the 
magnitude of the internal magnetic fleld, see Eq. ()5.1|) . 
we obtain the following condition when the separation 
into slow and fast variables is legitimate. The criterium 
can be formulated as Tq 3> t, which can be presented as 



^1 =sl 



(5.7) 



B. Applicability of the gaussian approximation 

Let us discuss the applicability of the gaussian approx- 
imation for the action in and h'^. The coefficients in 

front of terms hfl{t) h\\ (ti) . . . ft,|j'(t„) are obtained by tak- 
ing the rtth variational derivative of 5g^ {t, t) + 5g^ {t, t), 
or, equivalently, by iterating the Usadel equation n times. 
Since the typical frequencies of are uj ^ 1/r, the left- 
hand side of the equation is ~ /t, while the 
right-hand side is Since the only time scale 
here is r, all the coefficients of the expansion of the action 
in /i|j'(w) at w '-^ 1/t are of the same order: 

Q '^"^"^ f ■ ■ ■ d'^n 

X hfitoi) . . . hfiiJn) h^i-uji - ... - w„). (5.5) 



VI. MAGNETIC SUSCEPTIBILITY OF 
METALLIC PARTICLES OUT OF EQUILIBRIUM 

The LLG equation derived in this paper for a ferro- 
magnetic particle with finite bias between the leads can 
be applied to a number of experimental setups. More- 
over, the derivation of the equation can be generalized 
to spin-anisotropic contacts with leads or Hamiltonian of 
electron states in the particle. In this paper we apply 
the stochastic equation for spin distribution function to 
the analysis of the magnetic susceptibility at finite fre- 
quency. The susceptibility is the basic characteristic of 
magnetic systems, it can often be measured directly, and 
determines other measurable quantities. 

Below, we calculate the susceptibility of an ensemble 
of particles placed in constant magnetic field of an ar- 
bitrary strength and oscillating weak magnetic field, see 
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Fig. 1. We consider the oscillating magnetic field with 
its components in directions parallel and perpendicular 
to the constant magnetic field. 



A. Solution at zero noise power 

At Teff = when w{t) — 0, and at fixed direction 
of the field, B{t) = ezB{t), equation of motion (|4.4|) is 
easily integrated for an arbitrary time dependence B{t): 



The linear correction to the probability distribution can 
be cast in the form 

V{9, t)=[l + {6) e-'"* + e'"*] Vo{0) , 

(6.8) 

with Vo{9) defined by Eq. (|6.3p . The magnetic ac sus- 
ceptibility can be evaluated from Eq. (|6.8p as 



(6.9) 



X||(w,6) = 27r / u\i{9)Voie)cos9smede . 



if = ipo+ 2B{t') dt', 



tan — = tan — exp 
2 2 ^ 



hoT 



(6.1) 



(6.2) 



Here the direction of magnetic field corresponds to 9 = 0. 



B. Constant magnetic field 

At finite Tcff in constant magnetic field Bq the Fokker- 
Plank equation has a simple solution 



Vo{9) = 



sinh b 47r 



(6.3) 



where the strength of constant magnetic field is written 
in terms of the dimcnsionless parameter 



{2hoT)Bo 



(6.4) 



Substituting this probability function to Eq. (I4.15p . we 
obtain the classical Langevin expression for the magne- 
tization of a particle in a magnetic field 



M, = coth b - 



1 

6 ' 



^ My = 0. 



(6.5) 



This expression for the magnetization coincides with the 
magnetization of a metallic particle in thermal equilib- 
rium, provided that the temperature is replaced by the 
effective temperature T^s defined by Eq. (|3.30p . 
The differential dc susceptibility is equal to 



dc _ dM.jb) 
^11 db 



1 



sinh^ b 



(6.6) 



C. Longitudinal susceptibility 



We now consider the response of the magnetization to 
weak oscillations 13^ (t) of the external magnetic field with 
frequency lo in direction parallel to the fixed magnetic 
field Bq. We write the oscillatory component of the field 
in terms of the dimcnsionless field strength: 



bii e 



2hoB,{t) 
SiTcs 



(6.7) 



The equation for My (9) is obtained from Eq. (j4.10p with 
B, = Bn + B,{t): 

d'^UH cos9 - bsin^ 9 du\\ ^ , . 9^ „ 

+ — — f +zr!u|| =5sm2 0-2cos0, 

smO ad " 

(6.10) 

where we introduced the dimcnsionless frequency 

n^uj%, (6.11) 



and the time constant Tq is defined in Eq. (|4.9p . 

Note the symmetry of Eq. (|6.10p with respect to the 
simultaneous change b ~b and 9 ^ n — 9. Also, the 
normalization condition for the probability function re- 
quires that 



u\\iO)Voi9) sin9d9 = 0. 



(6.12) 



The latter holds if the boundary conditions Eq. (|4.14p 
are satisfied, which in the case of axial symmetry can be 
written as 



lim < sm 



dun(9) 
89 



0. 



(6.13) 



The differential equation (|6.10p with the boundary con- 
dition Eq. (j6.13p can be solved numerically and then the 
susceptibility is evaluated according to Eq. (|6.9p . The 
result is shown in Figs. [2] and [3l where the susceptibility 
is shown as a function of frequency w or magnetic field 6, 
respectively. We also consider various asymptotes for the 
ac susceptibility, obtained from the solution of Eq. (j6.10p . 

At zero constant magnetic field, 6 = 0, we find the 
exact solution of Eq. (|6.10p explicitly: 



cos 6* 



(6.14) 



This solution allows us to calculate the ac susceptibility 
in the form 



X„(^^,6 = 0) = l^ 



n/2 ■ 



(6.15) 



For 6^1 only cos 6* ^1/6 matter, and we can find a 
specific solution of the inhomogeneous equation: 



1 - 5(1 - cos( 
b~ifl/2 



6> 1. 



(6.16) 
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FIG. 2: (Color online). Plot of the real and imaginary parts 
of the susceptibility x\\ (^i as a function of the dimension- 
less frequency = ojTq. The oscillatory field at frequency uj 
is parallel to the constant magnetic field with strength 6. The 
real part of the susceptibility decreases monotonically from 
its dc value, Eq. (|6.6|l . as frequency increases, while the imag- 
inary part increases linearly at small Q, <^ 1, see Eq. (|6.20|l . 
and decreases at higher frequencies. 



The requirement of regularity at the opposite end can 
be replaced by the probability normalization condition, 
Eq. (|6.12p . which is satisfied by this solution. Substitut- 
ing this solution to Eq. (|6.9p . we obtain the strong field 
asymptote for the ac susceptibility 



1 



b{b-in/2) ■ 



(6.17) 



For ^ 1 and fl ^ b, we can neglect the derivatives 
in Eq. (|6.10p and find the solution in the form 



Me) 



b sin^ 6* - 2 cos ( 

in 



(6.18) 



This solution u\\{0) also satisfies Eq. (|6.12p . For the sus- 
ceptibility, Eq. (|6.9[) . we obtain 



XII (f^ ~* oo,6) 



2i / coth& 

^ n [ b 



1 

52 



(6.19) 



Finally, the low frequency limit can be also analyzed 
analytically. The real part of the susceptibility coincides 
with the differential susceptibility in dc magnetic field, 
Eq. (j6.6p . for the imaginary part to the first order in 
frequency we obtain, see Appendix, 



Imx||(r!,6) = r!/||(6). 



(6.20) 



The function /|| (b) has a complicated analytical form and 
is not presented here, but its plot is shown in Fig. [H 

In all considered four limiting cases, the asymptotic 
approximations hold regardless the order in which the 
limits are taken. Indeed, the asymptote of the expression 
for the susceptibility in the zero field, Eq. (|6.15p . has the 
asymptote at f2 ^ 00 consistent with Eq. (|6.19p at 6 = 0. 



FIG. 3: (Color online). Plot of the real and imaginary parts 
of the ac susceptibility x\\ b) at several values of the dimen- 
sionless frequency Q of the oscillating magnetic field along the 
constant magnetic field with strength b. In general, magnetic 
field suppresses both real and imaginary parts of the suscep- 
tibility. 



Similarly, the high frequency limit of Eq. (|6.17p coincides 
with the limit 6 — > cx) of Eq. (|6.19p . Both limits of weak 
and strong magnetic field of the imaginary part of the 
susceptibility at low frequencies, Eq. (|6.20p . coincide with 
the imaginary part of X|| b), calculated from Eq. (|6.15p 
and Eq. (|6.19p . respectively. 

In general, we make a conjecture that the ac suscepti- 
bility is given by the following expression: 



1 - in/Ti{b) ' 



(6.21) 



where functions Xn(^) ^^nd Tn{b) are real and describe 
the degeneracy points of the homogeneous differential 
equation Eq. (j6.10p with real iVL. This expansion is re- 
lated to the expansion of time-dependent Fokker-Plank 
equations in the spherical harmonics, analyzed in Ref. (H. 
In particular, xi>i(fe ^ 0) = 0{b) and vl{b ^ 0) = 
n{n + l) + 0{b). 

For practical purposes, we found from a numerical 
analysis that even the single-pole approximation. 



sinh^ b 



^ [l-*c.7]|(&)] \ (6.22) 



gives a very good estimate of the susceptibility for all 
values of lo and b. The analysis shows that the suscep- 
tibility, Eq. (|6.9p . obtained from a numerical solution of 
Eq. (|6.10p . is within a few per cent of the estimate given 
by Eq. (|6.22p . The characteristic time constant, 7j|(6), 
as a function of magnetic field b is chosen from the high 
frequency asymptote Eq. (j6.19p : 



r\\{b) = 



To 



sinh^ b — b^ 



2 sinh b b cosh b — sinh b 



(6.23) 



To evaluate the accuracy of the above approximation, 
Eq. (j6.22p . we consider the opposite limit of low frequen- 
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frib) 




FIG. 4: (Color online). Dependence on magnetic field b of 
the slope of the imaginary part of the linear in frequency 
susceptibility at low frequencies <C 1, calculated 

according to Eq. (|6.20|) . For comparison, we also plot function 
see Eq. mm . 



cies, n <JC 1, and compare the exact result for the imagi- 
nary part of the susceptibility, Eq. (|6.20p , with 



frib) 



1 



(sinh^ b~b^) 



2\2 



(6.24) 



2fo2 sinh b b cosh b - sinh b 



For visual comparison of functions /||(6) and /|'|''^''(6), we 
plot both functions in Fig. [4l where these curves are 
nearly indistinguishable. The difference between these 
two curves vanishes at 6 — > and b —^ oo, and has a 
maximal difference at & « 2, which constitutes only tiny 
fraction of /|| (6). 



D. Transverse susceptibility 

Next, we consider the response of the magnetization 
to weak oscillations B± (t) of the external magnetic field 
with frequency lo in direction perpendicular to the fixed 
magnetic field Bq. We write the oscillatory component 
of the field in the form: 



B±it) 



iJcff 



[b^{^ 



le. 



^)e— * + 61(e. 



(6.25) 

This field represents a circular polarization of an ac 
magnetic field in the {x, y) plane, perpendicular to 
the fixed magnetic field in the z-direction: B = 
{B ± cos ut; B± sin ujt; Bf)}. We look for the linear cor- 
rection to the probability distribution in the form 



x[l + b^u^{9) e^'^-*"* -f b\u\{e) 6"*'^+*'^*] 



(6.26) 



The equation for u±^{9) is obtained from the Fokker- 
Plank equation Eq. (|4.10p . linearized in the parame- 



ter b±: 
902 ■ 



cos — b s\T? 9 du I 



1 



smW 89 ' V"""" sin^9 

sin 0(2 + 2ihQTb + 6 cos 9) . 



Uj_ 



(6.27) 



Here the dimensionless frequency is a difference between 
the drive frequency uj and the precession frequency in 
external field Bq: 



nj_ = {uj- 2Bo) To = n- 2(/ior)&, 



(6.28) 



where Tq is defined in Eq. (|4.9p and the right equality is 
written in terms of dimensionless variables f2, Eq. (|6.11|) . 
and 6, Eq. ()6.4|) . Equation (|6.27p is symmetric with re- 
spect to the simultaneous change 9 tt — 9, 6 — > —6, 
i —i, il± — il_L ("parity"). The function 'P{(p,9,t) 
is single-valued at the poles 9 — and 9 = tt, only if 



u^{9 = 0) = 0, ujl{9 = tt) = 0. 



(6.29) 



The latter equations establish the boundary conditions 
for the differential equation (|6.27p . We also note that 
the normalization condition is satisfied for any function 
u^{9). 

We define the susceptibility in response to the ac mag- 
netic field, Eq. (|6.25p . as 

7r 

X±i^,b) = 2TT J u± {9)Vo (9) sin^ 9d9 . (6.30) 



This expression for the susceptibility can be used to cal- 
culate the magnetization of a particle 



Mit) 



n{ip, 9)V{ip, 9, t) sin 9 d9d(p 



(6.31) 



to the lowest order in the ac magnetic field. In particular. 



A'Ut) - Rc(x±6xe-*"*), My{t) = lm{x±b^e-''^'). 



..32) 



Solving numerically the differential equation (|6.27p 
with the corresponding boundary conditions, Eq. (|6.29p . 
we obtain the transverse susceptibility, Eq. (|6.30p . shown 
in Figs. [5] and Below we analyze several limiting cases. 

In zero fixed magnetic field, 6 = 0, we have the exact 
solution of Eq. ((07)l : 



u±{9) 



sin 9 



1 - inj_/2 



(6.33) 



This solution corresponds to the solution in the longitu- 
dinal case, rotated by 90°, cf. Eq. (pH)) . 

At w = and n± ^ -2bhoT, the solution of Eq. (jOTl) 
has a simple form U|| (9) — sin9 and corresponds to a tilt 
of the external field. The susceptibility due to such tilt 
is 



X±{^^0,b) 



62 



{bcothb- 1) . 



(6.34) 



13 
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FIG. 5: (Color online). Plot of the real and imaginary parts 
of the transverse susceptibility x±i^jb) as a function of the 
dimensionless frequency Q. Negative frequency corresponds 
to the opposite sense of the circular polarization of the ac 
magnetic field in a plane, perpendicular to the constant mag- 
netic field with strength b. The parameters of the three shown 
curves are chosen so that horb — 2. 



In strong fixed magnetic field, 6 ^ 1, we need to con- 
sider small angles 9 ^ ^/Vb, therefore, we can approxi- 
mate cos 9 1 in Eq. (|6.27p and obtain: 



b + 2 + 2ihoTb . 
b + 2 ~ iil± 

The susceptibility in the limit 6 3> 1 is given by 

1 + 2ihoT 



6(6(1 + 2ihoT) - ifl) 



(6.35) 



(6.36) 



At n± ^ 1, 6 we can disregard the terms in Eq. (|6.27p 
with derivatives. Moreover, the contribution to the sus- 
ceptibility, Eq. (|6.30p . from the vicinity of 6* = and 
9 = TT is suppressed as sin'^ 9. This observation allows us 
to write the solution in the form 



u±{9) 



smf 



2 + 2ihoTb + bcos9 



(6.37) 



Consequently, we obtain the following high frequency, 
n± ^ 1, asymptote for the susceptibility: 



X±{^,b) 



n - 2hoTb 



b coth 6—1 



(6.38) 

We can use the approximate expression for the suscep- 
tibility in response to the transverse oscillating magnetic 
field 

xrH-^^^^[i + 2.5or.(6)] 

X [1 + i{2Bo - uj)T^{b)r^ 
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FIG. 6: (Color online). Plot of the real and imaginary parts 
of the transverse susceptibility x±(f^jfe) as a function of the 
strength 6 of a constant magnetic field, shown for two values 
of frequency f2 and two values of the "damping factor" /lor. 
Negative values of 6 corresponds to the opposite sense of the 
circular polarization of the ac magnetic field in a plane, per- 
pendicular to the constant magnetic field. The real part of 
the susceptibility exhibits a strong non-monotonic behavior 
at weak magnetic fields. 



The corresponding characteristic time constant 71 (6) can 
be found for any 6 from the asymptotic behavior of 
X±(fi,6) at f}_L > 1,6: 



6 coth 6—1 
62 + 1 - 6coth6 



(6.39) 



VII. CONCLUSIONS 

We have studied the slow dynamics of magnetization 
in a small metallic particle (quantum dot), where the fer- 
romagnetism has arisen as a consequence of Stoner insta- 
bility. The particle is connected to non-magnetic electron 
reservoirs. A finite bias is applied between the reservoirs, 
thus bringing the whole electron system away from equi- 
librium. The exchange of electrons between the reser- 
voirs and the particle results in the Gilbert damping-^ of 
the magnetization dynamics and in a temperature- and 
bias-driven Brownian motion of the direction of the par- 
ticle magnetization. Analysis of magnetization dynam- 
ics and transport properties of ferromagnetic nanoparti- 
cles is commonly performed^i^i^iiiii within the stochas- 
tic Landau-Lifshitz-Gilbert (LLG) equation^'^, which is 
an analogue of the Langevin equation written for a unit 
three-dimensional vector. 

We derived the stochastic LLG equation from a mi- 
croscopic starting point and established conditions under 
which the description of the magnetization of a ferromag- 
netic metallic particle by this equation is applicable. We 
concluded that the applicability of the LLG equation for 
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a ferromagnetic particle is set by three independent crite- 
ria. (1) The contact resistance should be low compared to 
the resistance quantum, which is equivalent to Nch ^ 1- 
Otherwise the noise cannot be considered gaussian. Each 
channel can be viewed as an independent source of noise 
and only the contribution of many channels results in 
the gaussian noise by virtue of the central limit theorem 
for Nch ^1- (2) The system should not be too close 
to the Stoner instability: the mean-field value of the to- 
tal spin Sq ^ Nch- Otherwise, the fluctuations of the 
absolute value of the magnetization become of the or- 
der of the magnetization itself. (3) Sq ^ Toff/5i, where 
Teff ~ max{T, \ eV\} is the effective temperature of the 
system, which is the energy scale of the electronic dis- 
tribution function. Otherwise, the separation into slow 
(the direction of the magnetization) and fast (the elec- 
tron dynamics and the magnitude of the magnetization) 
degrees of freedom is not possible. 

Under the above conditions, the dynamics of the mag- 
netization is described in terms of the stochastic LLG 
equation with the power of Langevin forces determined 
by the effective temperature of the system. The effective 
temperature is the characteristic energy scale of the elec- 
tronic distribution function in the particle determined by 
a combination of the temperature and the bias voltage. 
In fact, for a considered here system with non-magnetic 
contacts between non-magnetic reservoirs and a ferro- 
magnetic particle the power of the Langevin forces is 
proportional to the low-frequency noise of total charge 
current through the particle. We further reduced the 
stochastic LLG equation to the Fokker-Planck equation 
for a unit vector, corresponding to the direction of the 
magnetization of the particle. The Fokker-Plank equa- 
tion can be used to describe time evolution of the distri- 
bution of the direction of magnetization in the presence 
of time-dependent magnetic fields and voltage bias. 

As an example of application of the Fokker-Plank 
equation for the magnetization, we have calculated the 
frequency-dependent magnetic susceptibility of the par- 
ticle in a constant external magnetic field (i. e., linear 
response of the magnetization to a small periodic mod- 
ulation of the field, relevant for ferromagnetic resonance 
measurements). We have not been able to obtain an ex- 
plicit analytical expression for the susceptibility at ar- 
bitrary value of the applied external field and frequency; 
however, analysis of different limiting cases has lead us to 
a simple analytical expression which gives a good agree- 
ment with the numerical solution of the Fokker-Planck 
equation. 
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APPENDIX A: LONGITUDINAL 
SUSCEPTIBILITY AT LOW FREQUENCIES 

We find the linear in frequency 17 1 correction to the 
dc susceptibility. For this purpose, we look for a solution 
to Eq. (pl^ in the form 



\{e) = u^\e) + u'^>{e) 



(1)/ 



(Al) 



where u^«'\9) is the solution of Eq. (|6.10p at 17 = and 



u^^\0) cx 12. We choose 



,(0) 



{0) = \- coth b -+ cos e, 



(A2) 



-II h 

since this form of u"^^ (9) preserves the normalization con- 
dition (|6.12|) . This function can be found directly as a 
solution of Eq. (|6.10p with 17 = or as a variational 
derivative of function Vo{0), defined in Eq. (|6.3p . with 
respect to b. 

The linear in 17 correction u^^\6) is the solution to the 
differential equation 



(1) 



96*2 sin 61 09 

From this equation, we can easily find 



= -tnuf\0). (A3) 



du\^\e) 

89 



in 

6sin( 



-b cos 6 



coth b — COS 9 



sinh6 



(A4) 



We notice that the solution to the latter equation will 
automatically satisfy the boundary conditions, given by 
Eq. ()6.13p . Integrating Eq. (jA4p once again, we obtain 
the following expression for function uiP{9): 



u^\9)^C{b) 



coth b — cos 9' — 



-b cos 8' 



sinh b 



d9' (A5) 



smf 



Here the integration constant C{b) has to be chosen to 
satisfy the normalization condition, Eq. (|6.12p , which re- 
sults in complicated expression for the final form of the 
function u^^^'"' (9). 

To obtain function /||(6), introduced in Eq. (|6.20p . we 
have to perform the final integration 



/||(^) = 



2n 

"it" 



u\^\9)roi9) sine cos 



(A6) 



The result of integration is shown in Fig. |4l 



15 



1 W. F. Brown, Jr., Phys. Rev. 130, 1677 (1963). 
^ L. Landau and E. Lifshitz, Phys. Z. Sowietunion 8, 153 
(1935). 

3 T. Gilbert, Phys. Rev. 100, 1243 (1955). 

* J. L. Garcia-Palacios and F. J. Lazaro, Phys. Rev. B 58, 

14937 (1998). 
5 K. D. IJsadcl, Phys. Rev. B 73, 212405 (2006). 
® J. Foros, A. Brataas, G. E. W. Bauer, and Y. Tserkovnyak, 

Phys. Rev. B 75, 092405 (2007). 

S. I. Denisov, K. Sakmann, P. Talkner, and Hanggi, Phys. 
Rev. B 75, 184432 (2007). 
® A. Rebei and M. Simoniato, Phys. Rev. B 71, 174415 
(2005). 

^ H. Katsura, A. V. Balatsky, Z. Nussinov, and N. Nagaosa, 

Phys. Rev. B 73, 212501 (2006). 
^° A. S. Niinez and R. A. Duine, Phys. Rev. B 77, 054401 
(2008). 

A. L. Chudnovskiy, S. Swiebodzinski, and A. Kamenev, 

Phys. Rev. Lett. 101, 066601 (2008). 

J. Foros, A. Brataas, Y. Tserkovnyak, and G. E. W. Bauer, 



Phys. Rev. Lett. 95, 016601 (2005). 

R. A. Duine, A. S. Niifiez, J. Sinova, and A. H. MacDonald, 
Phys. Rev. B 75, 214420 (2007). 

X. Waintal and P. W. Brouwer, Phys. Rev. Lett. 91, 

247201 (2003). 

I. L. Aleiner, P. W. Brouwer, and L. I. Glaznian, Phys. 
Reports 358, 309 (2002). 

C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 
I. L. Kurland, I. L. Aleiner, and B. L. Altshuler, Phys. 
Rev. B 62, 14886 (2000). 

Y. Tserkovnyak, A. Brataas, and G. E. W. Bauer, Phys. 

Rev. Lett. 88, 117601 (2002). 

S. Tarnaru et ai, ,J. Appl. Phys. 91, 8034 (2002). 

J. C. Sankey et ai, Phys. Rev. Lett. 96, 227601 (2006). 

Y. Ahmadian, G. Catelani, and I. L. Aleiner, Phys. Rev. 

B 72, 245315 (2005). 

K. D. Usadel, Phys. Rev. Lett. 25, 507 (1970). 

G. Catelani and M. G. Vavilov, Phys. Rev. B 76, 201303 

(2007). 



